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Abstract 



The equation of state of pure QCD, obtained from lattice QCD, is discussed 
for temperatures ranging from 0.9T C to 4T C , as well as results on screening 
masses, the chiral condensate, and the pion decay constant close to the de- 
confinement phase transition in the confined phase of QCD. The equation of 
■ state differs significantly from that of a free gas. There is little evidence of a 

temperature dependence in the chiral condensate or the meson properties, but 
perhaps some for the nucleon screening mass. Above the phase transition one 
sees non-perturbative effects, even though hadron correlators show the exis- 
tence of deconfined quarks. 

o 
in 
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1 Introduction 

The study of hadronic matter at high temperature and / or density is receiving in- 



creasing attention, both experimentally and theoretically. It is important to gain 
a good understanding of the predictions of QCD in this regime as soon as possi- 
ble. There are many approaches to this task, some of which are presented elsewhere 
in these proceedings. Each is advantageous for some aspect of the problem. One 
needs to combine them all to understand what QCD predicts over the entire range of 
temperatures and densities relevant for heavy ion collisions. 

Recent results from one of the major theoretical players, the lattice regulation 
of QCD at finite temperature, will be presented below. An introduction to lattice 
QCD may be found, for example, in [0, 0. 

When one is close to the phase transition, many of the traditional approaches to 
QCD run into problems, as both low temperature and high temperature expansions 
need to be extrapolated into regions far from the regions in which they are known to 
be secure. The lattice may well be the only satisfactory approach in this region. 

The equation of state, both for nuclear matter just below the deconfinement 
transition, and for the quark-gluon plasma above the phase transition, is required 
information for, amongst others, hydrodynamic models attempting to describe the 
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evolution of the quark gluon plasma. Since the couplings and masses are not small 
below 2T C the lattice approach is the obvious method to use. 

Despite extensive studies of the QCD phase transition [Q, |3| little is known with 
certainty about the changes in the excitation spectrum. In particular, the properties 
of the (quasi)-particle spectrum in the confined phase close to T c require a more 
thorough understanding, as the temperature dependence of hadron masses and other 
hadronic parameters will lead to observable consequences in current, and especially 
forthcoming, heavy ion collisions. 

The temperature dependence of various hadronic properties has been addressed 
using a number of different approaches. These calculations yield different predictions 
for the behaviour of some quantities |4| . Since these approaches have limited applica- 
bility at high temperatures, one may hope that a lattice calculation of the variation 
with temperature of any of these quantities below the phase transition may shed light 
on the discrepancies. 

Another important quantity is the ratio of the chromo electric to chromo mag- 
netic screening length. This is used, for example, in estimates of the colour trans- 
parency of the plasma [||. One promising approach is to calculate it directly from 
measurements of the gluon self energy on the lattice. 

Recent results on the equation of state for pure gluon SU(3) theory and for 
QCD with two light flavours will be presented in section |2|. The results from a study 
of hadronic properties in quenched QCD at temperatures between 0.75T C and 0.92T C 
are presented in section Section f| contains the results pertaining to the chromo 
electric and magnetic screening masses. 



2 Equation of State 

Recent work || has enabled the calculation of the energy density and pressure on the 
lattice using an entirely non-perturbative technique for the first time. This requires 
detailed a knowledge of dg~ 2 /da, ie. one needs to know the QCD beta function /3f 
relating changes in physical quantities with changes in the coupling. This has been 
calculated by the TARO collaboration J7]. 

Given the beta function one calculates the interaction measure, the difference 
between the energy density and thrice the pressure, as follows: 



e — 3p 

rp4 



T— (JL 

dT \T 4 
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The pressure is given by 
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(2.2) 



The zero temperature plaquette is given by Po, the spatial and temporal plaquettes 
at finite temperature are given by P a and P T respectively. Both the energy density 
and the pressure are normalised to zero at zero temperature. 
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Figure 1: The interaction measure (e — 
3P)/T 4 against temperature for three 
different values of the temporal extent N T 
in pure SU(3) gauge theory. 



Figure 2: The energy density compared 
to the energy density in the Stefan- 
Boltzmann limit on a lattice of the same 
size against the temperature in units of 
the critical temperature. 



The results obtained for the interaction measure, divided by T 4 , are shown in 
figure [[} The peak in the interaction measure shortly after the first order deconfine- 
ment phase transition shows clearly the strong deviation from ideal gas behaviour, 
e = 3P. The energy density rises sharply, whilst the pressure lags behind until one 
has reached temperatures around 3T C . This is evident again in figure |3] where the 
energy density and pressure, scaled with T 4 are plotted together. 

Note that the interaction measure is related to the value of the gluon condensate 
at finite temperature f| || via 

(G 2 ) T =(0|G 2 |0)-(e-3P) T , (2.3) 

where {G 2 )t and (0|G 2 |0) are the condensates at zero and finite temperature respec- 
tively. The gluon field strength G 2 is given by G 2 = -[/V(2# 3 )]GTG%. Taking 
(0|G 2 |0) = 0.0135GeV 4 f| one finds that (G 2 ) T drops to zero immediately after the 
phase transition, and then goes negative according to the above equation. 

The energy density for different values of the temporal extent is plotted in 
figure |2], compared to the Stefan-Boltzmann limit on this size lattice. If one takes the 
continuum limit one obtains a value around 83% of the Stefan-Boltzmann value for 
temperatures above 2T C up to at least 3T C . It is clear that even at these temperatures 
the interactions and masses of the excitations cannot be neglected. A similar rapid 
rise in the energy density, and slower rise in the pressure, is seen in the case of two 
flavour QCD JTCj . 

It is clear that, if there is a jump in the energy density, the pressure will lag 
behind in a thermodynamic system, from the relation P/T = Pq/T + / dTe/T 2 . 
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Figure 3: The energy density and pres- 
sure in units of T 4 on a lattice of size 
32 3 x 6. Also shown is the continuum 
Stefan-Boltzmann limit (lower solid line) 
and the Stefan-Boltzmann limit for this 
size lattice (upper solid line). 



Figure 4: The speed of sound squared, 
Cg = dP/de function of temperature 
on lattices of size 16 3 x 4 and 32 3 x 6 as a 
function of temperature. 



The deviation from the ideal gas behaviour may be understood in terms of plasma 



consitituents that are massive, and have a strong interaction [] I 



The speed of sound in a medium is another important parameter. This is shown 
in figure f| for the high temperature phase of pure gauge theory. The speed of sound 
also approaches the ideal gas value (l/v3), but is only at around 80% of it even at 
2T r . 



3 Hadronic Properties 
3.1 Chiral Sector 

The quark condensate, the quantity from which other approaches obtain much of the 
effect of temperature on the hadrons, is the logical place to start. It is known to 
change drastically at the deconfinement transition in quenched QCD. 

The chiral condensate is shown in figure [|(a) as a function of bare quark mass 
at T = 0.92T C . It may be calculated directly from the trace of the fermion matrix 
((■0'0)se) and by using the pion and sigma propagators, which corrects for the linear 



dependence on the quark mass [pl|. One can see quite clearly that both methods 
extrapolate to the same value at zero quark mass. The same holds true for lower 
temperatures. 

The temperature dependence of the chiral condensate is shown in figure H(b), 
where the ratio of the finite temperature chiral condensate to that at T = is plotted. 
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Figure 5: Figure (a) shows the chiral condensate as a function of the bare quark 
mass, at (3 = 6.00. The upper curve is obtained using the value taken directly from the 
trace of the fermion matrix, the lower curve after subtracting the linear dependence 
on the quark mass. The results have been taken from [12|, 13]. Figure (b) shows 



the finite temperature chiral condensate normalised to its value at zero temperature. 
The circles represent quenched data, extrapolated to zero quark mass from fl2[] . The 
other figures represent results for various numbers of flavours, but not extrapolated 
to zero quark mass. 



These values have been extrapolated to zero quark mass, which is the reason why the 
error bars are so large. If one does not extrapolate, the ratio remains the same, but 
with errors a factor of ten smaller. For comparison results for full QCD, including 
dynamical fermions, have been included. Since none of these results are at the physical 
quark mass, the effect of temperature is probably underestimated. It is clear that 
the chiral condensate does not change significantly until one is extremely close to the 
critical temperature. 

The pion mass displays the behaviour expected of a Goldstone particle, = 
A n m q , at all temperatures. There is also no sign of a dependence on temperature in 
the pion mass, or in the slope A v . 

Another quantity of considerable interest for understanding the temperature de- 
pendence of the chiral sector of QCD is the pion decay constant. It can be determined 



directly from the relevant matrix element on the lattice |i| . At zero temperature f n 



is related to the chiral condensate and the pion mass through the Gell-Mann-Oakes- 
Renner relation, which is expected to be valid as long as chiral symmetry remains 
spontaneously broken. 

m lfl = "lq<W(mq=0) • (3-4) 

The temperature independence of the amplitude, A n , and the chiral condensate 
up to T = 0.92T C indicates that the pion decay constant will also be temperature 
independent. A direct calculation shows that it is, and that the GMOR relation 
holds, for temperatures up to at least T = 0.92T C . 
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3.2 Nucleon and meson masses 



One obtains hadron screening masses by fitting the exponential decrease of the corre- 
lator Cr(z) to, for example, an hyperbolic cosine. One obtains either local masses, if 
only two points are fitted to, or an estimate of the lowest state if one fits to the long 
distance part of the correlator. The lowest mass which can be extracted from the fit, 
.Eh, is related to the screening mass for fermionic states via 

E^ = m 2 n + k sin 2 (n/N T ) , (3.5) 

due to the contribution from the non-vanishing Matsubara energy, po = vrT. Note 
that for bosons = m H . 

The local rho screening masses and the value obtained from a fit to the full 
propagator at T = 0.92T C and at zero temperature are shown in figure ^|(a). It 
is clear that the rho screening mass does not show any significant dependence on 
temperature. 

The nucleon has been examined using wall sources at a quark mass of 0.01 in 
lattice units. The local screening masses for T = 0.92T C , and the screening mass 
obtained from a full fit, are shown in figure ||(b) along with the result obtained at 
zero temperature. The increase in temperature clearly has a dramatic effect on the 
nucleon. However, most of this can be understood purely in terms of the fact that 
the lowest momentum of the nucleon is not zero, but irT, as discussed above. 

The screening mass at T = 0.92T C , extracted using eqn. |3]5| with the assumption 
that k = 1, indicates that the mass rises slightly: m n (0.92T c ) = (l.l±0.03)m n (T = 0). 
We note, however, that there is also the possibility of a modification of the energy 
dispersion relation at finite temperature, which may lead to a deviation of k from unity 
|2"0fl . The error given above is from the statistical error alone, as we cannot determine 



the systematic error which may be introduced by our assuming that eq. |3.5| with k — 1 
is applicable. 

Finally, one would like to know the properties of correlators with hadronic 
quantum numbers in the deconfined phase. These can be studied using correlators in 
the spatial direction, as above. One finds that the nucleon then approches thrice the 
lowest Matsubara frequency, and all mesons other than the pion approach twice the 
lowest Mastubara frequency (see [H and refernces therein). This is what one expects 



if the quarks become deconfined. 

The pion, though, acquires a screening mass half as large as expected for a free 
meson. As it is the particle associated with chiral symmetry, one may expect it to have 
a special role even above the phase transition [^|. In order to understand whether 
the pion is indeed deconfined as well, with a different interaction between the two 
quarks, the meson correlators were examined using different techniques ||23|| . In one 



of these, temporal correlators generated with wall sources were examined. Temporal 
correlators show the true mass, and not the screening mass. Wall sources are used in 
order to project out only the contribution from the lowest state. The results for the 
local masses are shown in Figure [7|, along with results for the effective quark mass. 

One sees here that both the rho and the pion masses tend towards twice the 
effective quark mass at high temperatures. This, combined with other investigations 
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Figure 6: Figure (a) shows the local rho screening mass at T = 0.92T C , and 
m q = 0.01 in lattice units, from wall sources. The solid and dotted lines show the 



fitted value and errors, the long dashed lines the zero temperature values from [ plf . 
Figure (b) shows the local nucleon screening mass at T = 0.92T C , and m q = 0.01 in 
lattice units, from wall sources. The solid and dotted lines show the fitted value and 



errors, the long dashed lines the zero temperature values from [21 



Figure 7: Local masses m H wa from correlators constructed with wall sources in 
the PS (circles) and V (squares) channels in the high temperature phase. Also shown 
is 2rriq S (triangles) from [Q. Lines have been drawn to guide the eye. 
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on the lattice, indicate that the hadrons (including the pions) become deconfined 
above T c , with rather large channel dependent residual interactions. 

One does, in all cases, see the factor of two difference between the pion and 
rho masses or screening masses. There is clearly a difference between a quark-anti- 
quark pair carrying pion quantum numbers and one carrying rho quantum numbers. 
Explanations for the relative lightness of the pion have been put forward based on 
spin-spin interactions |25| and residual gluon condensates p6 . 



4 Electric and Magnetic Screening Masses 

The chromo electric screening mass can be obtained by examining, amongst others, 
the correlation between Polyakov loops or the pole of the gluon correlator. At high 
temperature in perturbation theory at next to leading order it has been shown that 



the two approaches yield the same result |27, |28 



The chromo magnetic mass cannot be determined in perturbation theory, due 
to infra-red divergences. Preliminary results from calculations in progress 



m 



pure gauge theory, indicate that the magnetic and electric screening masses have the 
same order of magnitude above the phase transition, and for temperatures as high 
as 2T C . This result is consistent with the results obtained in another study in SU(2) 
pure gauge theory p3| . 

This can be compared with the results from the potential, and the spatial 
string tension |3T] which indicate that the coupling remains large well above the 



phase transition, with g(T = 2T C ) « 2. 



5 Conclusions 

The chiral condensate does not show any significant temperature dependence up to 
T = 0.92T C . In view of this the observation that there is no sign of a temperature 
dependence in the pion decay constant, pion screening mass, sigma screening mass or 
the rho screening mass may not come as a surprise. There is some, albeit inconclusive, 
indication of a temperature dependence for the screening mass of the nucleon. 

The details of the temperature dependence may change in the case of QCD 
with two flavours, where the transition is expected to be second order rather than 
first order as in the the pure SU(3) gauge theory, and with quark masses close to 
the physical quark mass. However, these differences will probably not show up for 
temperatures less than 0.9T C . 

The equation of state for pure gauge theory is now reasonably well understood. 
One finds clear evidence for a departure from the ideal gas at below 2T C , at which 
point the energy density reaches only 83% of the Stefan-Boltzmann limit. Thus it is 
clear that there are strong interactions, or equivalently a large contribution from non- 
perturbative effects, above the phase transition. This is supported by measurements 
of the spatial string tension, which indicate a coupling g w 2, results for the chromo 
electric and magnetic screening masses, as well as the effective quark mass in the 
Landau gauge. 
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All states with hadronic quantum numbers appear to be deconfined above the 
phase transition. However, the coupling between the quark and anti-quark depends 
on the channel considered, leading to pion / sigma masses of around half the rho 
mass. 
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